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ABSTRACT 

The violent giant flares of magnetars excite QPOs which persist for hundreds of 
seconds, as seen in the X-ray tail following the initial burst. Recent studies, based 
on single-fluid barotropic magnetar models, have suggested that the lower- frequency 
QPOs correspond to magneto-elastic oscillations of the star. The higher frequencies, 
however — in particular the strong 625 Hz peak — have proved harder to explain, 
except as high mode multipoles. In this work we study the time evolutions of non- 
axisymmetric oscillations of two- fluid Newtonian magnetars with no crust. We consider 
models with superfluid neutrons and normal protons, and poloidal and toroidal back- 
ground field configurations. We show that multi-fluid physics (composition-gradient 
stratification, entrainment) tends to increase Alfven mode frequencies significantly 
from their values in a single-fluid barotropic model. The higher-frequency magnetar 
QPOs may then be naturally interpreted as Alfven oscillations of the multi-fluid stel- 
lar core. The lower-frequency QPOs are less easily explained within our purely fluid 
core model, but we discuss the possibility that these are crustal modes. 
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1 INTRODUCTION 

The quasi-periodic oscillations (QPOs) of magnetars provide 
a tantalising possibility of probing the interior physics of 
neutron stars. Detected in the aftermath of giant flares, these 
■ QPOs are thought to relate directly to oscillation modes of 
$H ' the underlying star, perturbed by the hugely energetic flare. 
The theoretical challenge is to provide sufficiently good mod- 
els of magnetar oscillations that the observed QPOs may be 
convincingly linked to particular modes; this could provide 
potential constraints on the stellar equation of state and 
magnetic field configuration. 

Magnetars are thought to be slowly-rotating and 
magnetically-powered neutron stars (NSs), with dipole sur- 
face values reaching 10^^ G — an exceptional value, 
even for neutron stars. Three of these objects have 
been seen to pr oduce enormously energetic 'giant flares': 
SGR 0526-66 (iBarat et al ' 1983), SG R 1900-H4 and 



SGR 1806-20 dlsrael et al.ll2005. : Strohmaver fc Wattsll2005l : 
I Watts fc Strohmavei ]2006l ) . In the latter two cases, QPOs 



were detected in the flare's decaying X-ray tail, between 100 
and 400 seconds after the initial burst. A number of frequen- 
cies were identified within each tail, of varying duration and 
intensity. 
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Initially these QPOs were identified with crustal shear 
modes, and there was hope of using these to constrai n 
the equation of state (jPunca 3 Il998l : llsrael et al.1 \2004 ). 
Since these optimistic early studies, the picture has become 
more complex, due to the effect of the magnetic field: the 
crust and core are co upled, with r nodes becoming magneto- 
elastic in character iLevinl l2006l : iGlampedakis et aLlTiooel : 
iGabler et al.ll201^ . Attempts to model observ ed QPOs as 
global modes were dealt a blow bv lLevinl()2007l l. who used a 
toy model to show that axisymmetric oscillations of a mag- 
netic star could form a continuum of frequencies, with dis- 
crete crustal modes quickly damped and QPOs representing 
the edges of the Alfven continuum. Numerica l simulations 
supported this conclusion in the axial case (jSotani et al.l 
2008), but found that polar oscill ations were discrete modes , 
with no evidence for a continuum jSotani fc Kokkota J2OO9I ). 
These two studies considered axisymmetric oscillations on a 
purely poloidal background field. 

There are now indications that situations with coupling 
between axial and polar sectors may destroy the contin- 
uum. Two studies of non-axisymmetric magnetar oscilla- 
tions — for purely toroidal and purely poloidal background 
fields — found only discrete modes (iLander et al.l [2OI0I : 
iLander fc Jonesll201lh . These seemed at odds with the re- 
sults for axisymmetric studies; the key difference may be 
that non-axisymmetric Alfven modes are not purely axial 
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or purely polai0. For axisymmetric modes, a purely poloidal 
background field has decoupled axial and polar modes, but a 
mixed poloidal-toroidal field couples the two sectors. These 
coupled axial-polar oscillations were also recen tly found to 
be discrete modes l|Colaiuda fc Kokkotas|[2012l V 

Despite the great increase in sophistication of magnetar 
QPO modelling, a major missing ingredient in all these stud- 
ies has been the multi-fluid, non-barotropic nature of neu- 
tron star matter. A large fraction of mature neutron stars 
(including magnetars) is comp osed of superfluid neutrons 
and superconducting protons l|Bavm et al.l Il969l : IHo et all 
and the resulting oscillation spectrum will surely be 
different from that of a barotropic star. Furthermore, there 
are indications t hat no stable magnetic equilibria exist in 
barotropic stars (|Reisenegged I2OO9I : iLander fc Jon3 l2012l l 
— in which case it is clearly undesirable to use them as 
background models when studying oscillations. The rea- 
son why many magnetar QPO studies have not encoun- 
tered these instabilities so far is that they specialise to ax- 
isymmetric oscillations on poloidal-field backgrounds; the 
strongest instabilities appear for non-axisymmetric pertur- 
bations (Markev fc Ta ylor 1973; Wright 1973). 

In this work we try to advance the modelling of magne- 
tar QPOs, by studying the oscillations of a stratified multi- 
fluid star. We consider poloidal and toroidal magnetic field 
geometries. Two main simplifications we make are that the 
whole star is fiuid (with no crust), and that the neutrons 
are superfluid but the protons normal; the validity of these 
assumptions is discussed in the following sections. 

We begin by reviewing in Section [2] the main formal- 
ism for studying the dynamics of a magnetised two-fluid 
star. First we determine the background conflgurations and 
then we derive the linearised dynamical equations. In Sec- 
tion |3] we present our results, where we explore the differ- 
ences from barotropic single-fluid magnetar models. Finally 
in Section U we discuss our results in the context of magne- 
tar QPOs and potential future improvements. 



2 EQUATION OF MOTION 

The basic matter constituents of a neutron star are neutrons, 
protons and electrons (npe), although more exotic particle 
species may exist in an inner core region. In this work we 
assume a npe-composition in the entire volume of the star 
and neglect the inner core. 

From a dynamical point of view, we may assume that 
the core's protons and electrons can be considered as a sin- 
gle co-moving fluid, as the electromagnetic interaction locks 
them on timescales much smaller than the dynamical oscil- 
lation periods. The dynamics of this system can be therefore 
described by a two-fluid model comprising a gas of super- 
fluid neutrons and a neutral mixture of protons and electrons 
that we call for simplicity "protons". Neutron and proton 
quantities will be indicated with a subscript roman n and p, 
respectively. 

The equations to study the dynamics of a superfluid 
and magnetised star are, in the ideal MHD approximation, 

^ Instead, they are 'axial-led' o r 'polar-led', using the terminology 
of lLockitch fc Friedmaiil lll999l) . 



the mass and momentum conservation equations of each 
fluid constituent, the induction equation for the magnetic 
field, and the Poisson equa tion for the gravitational poten- 
tial. These equations read (iGlampedakis et al.ll20lil '): 

dtp^ + V • {pa) = , (1) 

F 

{dt + «xV) («" + e^wy^) + V{^ + fl^) + exu>f Vu^ = — , 

Px 

(2) 

dtB=V X {vpX B), (3) 
V^$ = 47rGp, (4) 

where the labels x and y (with x 7^ y) denote the fluid 
component n and p. The quantities px, Px and are, re- 
spectively, the mass density, chemical potential and velocity 
of each fluid constituent, while the relative velocity is de- 
noted by t^xy = Dx ~ I'y • The gravitational potential and the 
magnetic fleld are described by $ and B, respectively, and 
p = Pn -I- Pp is the total mass density. 

The parameter Ex accounts for the entrainment between 
nucleons, which is a non dissipative effect that in neutron 
stars is due to the strong interaction. The main effect of en- 
trainment is to couple the nucleon motion by inducing a rela- 
tive dragging between neutrons and protons. As a result the 
conjugate momentum of each component is not aligned with 
its velocity, see equation ([2]). The entrainment parameter 
can be also written in terms of the nucleo n's effective mass 
m* by using the relation £x = 1 — m*/mx l|Prix fc RieutordI 
2002). 

The vector field in equation ([2} represents the force 
density that acts on the x fluid component. The magnetic 
interaction and the mutual friction are the forces expected in 
a magnetised superfluid star with no crust. For normal (not 
superconducting) protons the interaction with the magnetic 
field is given by the Lorentz force: 

Fl = ^{VxB)xB, (5) 
47r 

where we have used Ampere's law to replace the charge 
density current. The mutual friction is instead a dissipa- 
tive force mediated by superfiuid vortices and operates in 
rotating stars on both the fiuid components. In this work 
we are interested in the oscillation spectrum of magnetars, 
which are slowly rotating objects that can be well described 
by non-rotating stellar models. We can therefore neglect the 
effects of rotation and mutual friction. More details on the 
impact of mutual friction on the o scillation spectrum in un- 
magne tised stars can be found in IPassamonti fc AnderssonI 
(|201ll ). 

Under these assumptions, the force density in equa- 
tion ((2]) is given by Fp = Fl and Fn = 0. 



2.1 Equation of State 

The equation of state (EoS) can be described by an energy 
functional 

£ = £ {pn,Pp,wlp) , (6) 
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that ensures Galilean invariance. The chemical potential /i 
and the entrainment parameter £x are then defined by 

9px 



(7) 



2px 



Py,W 

d£ 



(8) 



If the relative velocity between the two fluids is small, a 
typical situation in most astrophysical systems, equation ((6]) 
can be expanded in series 



£ ^ So {pn, Pp) + ao {pn, Pp) + O {wtp) 



(9) 



and the bulk EoS £o and the entrainment parameter qq 
can be independently specified at Wnp = 0. This approxi- 
mation is certainly valid in our star's model as a non zero 
relative velocity appears at first perturbation order. From 
equation (|8]l it follows that the entrainment parameter Ex is 
given by pxEx = 2ao. 

We consider an analytical Eo S, which is a two-fluid 
analogue of a p olytropic model ("Prix & Rieutord 2003: 
Andersson et al.ll2002 : .Passamonti et al.i,200a : .Lander et alj 
2012D : 



£o = kn Pn" + kp pp^ 



(10) 



where the coefficients fcx are constants and 7x is related 
to the adiabatic index A'x by the standard definition 7x = 
1 -I- 1/A^x- Despite its simplicity, this EoS allows us to con- 
struct neutron star models with composition gradients when 
A'n 7^ A'^p. Stratification is a relevant property for the oscil- 
lation spectrum and may be crucial for the stability of the 
background magnetic configuration. 



2.2 Background models: magnetic two-fluid 
stellar equilibria 

We begin by modelling a neutron star as an axisymmetric 
body in Newtonian gravity. For this system it is natural to 
use cylindrical polar coordinates {■uj,4>,z), where the z-axis 
is aligned with the symmetry axis of the star's magnetic 
field. 

We assume the whole interior of the magnetar is multi- 
fluid, with the neutrons superfluid, but the protons normal 
rather than superconducting. Although this is done for sim- 
plicity, it may be realised in magnetars if their internal field 
strength reaches ~ 10^^ G; at this 'upper critical fi eld' su- 
perconductivity is broken (|Glampedakis et al.ll201l] ). We do 
not account for the presence of the neutron star crust in this 
study. 

Since a full description of the equilibrium equa tions and 
their solution are presented in lLander et al.l (|2012l ). we pro- 
vide only a brief summary of the important details here. We 
neglect the inertia of the electrons and absord their chemical 
potential into that of the protons, so our equations describe 
a two-fluid system. From equation ([2]), we see that the sep- 
arate Euler equations governing the two fluids are: 



V(/in + $) =0, 



V (Ap + 



(V X B)xB 

47rpp 



(11) 
(12) 



for the gravitational potential, from which we see that the 
behaviour of the two fluids is coupled, despite the neutrons 
being a superfluid. 

The equilibrium magnetic field of an axisymmetric 
barotropic star is governed by a single equation in terms 
of the magnetic streamfunction u, known as the Grad- 
Shafranov equation. Although our stellar models are multi- 
fluid and stratifled, our chosen equation of state still allows 
us to derive a variant of the usual Grad-Shafranov equation, 
by replacing the total density p with the proton-fluid density 
pp- 

a' I d _^ ^ 2dM d/ 



vj dzu dz 



du 



where M = M{u) is a scalar function related to the Lorentz 
force and / — f{u) dictates the structure of the toroidal 
component in a mixed-fleld configuration. This equation de- 
scribes equilibria with purely poloidal and mixed poloidal- 
toroidal fields. We defer mixed fields to future work, and 
therefore set f{u) = to produce pure-poloidal field mod- 
els. Pure-toroidal field models are not governed by a separate 
equation but just involve an extra term in the proton-Euler 
equation. 

We employ an iterative numerical scheme to self- 
consistently solve the equations and H10p - (|13p in integral 
form; this produces our background equilibria. 

2.3 Perturbation Equations 

The perturbation equations for studying magnetised super- 
fluid stars can be derived by linearising equations H])-© 
and using as dynamical variables the mass density, veloc- 
ity and magnetic fleld perturbations. However, we prefer to 
define a new set of perturbation variables which both re- 
flects the dynamical degrees of freedom of a two-fluid sys- 
tem and simplifies the implementation of the boundary con- 
ditions. A pulsating two-fiuid star behaves dynamically as a 
coupled harmonic oscillator in which we may discern a co- 
and counter-moving relative motion between the two con- 
stituents. In general, these two degrees of freedom are cou- 
pled, but in some particular cases they can be completely 
decoupled — for instance in a non-stratified stellar model 

with no magnet ic field. 

Following iPassamonti et ah I (|2009l ). we use equa- 
tions ([l])-(l2]) to derive the following system of dynamical 
equations for non-rotating stars: 

dtf = -VSP+—5p-pVS$ + Pp5{Fl/Pp), (14) 
P 

dtD = (1 - xp) Pp [-S/S13 + 5 (Fl/pp)] , (15) 

dtSp = -V-/, (16) 

dtSxp = -V D-f Vxp. (17) 

where we have defined 7^ = 1 — — Cp and the following 
fluid perturbation variables: 



/ = Pn<5Wn + PpSVp , 

5p — Spn + Spp , 

D — Pp (1 — Xp) Swpn , 

^Xp = pSxp. 



(18) 
(19) 
(20) 
(21) 



As in the single-fluid case we have Poisson's equation Q The quantity SP is the total pressure perturbation that for 
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a co-rotating background reads 

V5P = 5 (pnV/in + PpV/ip) 



(22) 



while 5/3 = Sjlp — S/lu describes the deviation from beta equi- 
librium induced by oscillations and 5$ is the gravitational 
potential perturbation. 

To simplify the boundary conditions we choose a "flux" 
variable also for the magnetic field perturbation, i.e. we de- 
fine b = PpSB. The induction equation then reads 

dtb = XpV X (/x B) + V X (Dx B) - Xp'^ x (fx B) 

P 



x{Dx B) 



(23) 



while the perturbation of the Lorentz force in equations (|14p 
and (|15|) assumes the following form: 

&Fl = {V X B) X — + {V X b) X — - \ (Vpp xb)xB. 
Pp Pp Pp 

(24) 

It is worth noticing that the Alfven 

vel ocity in a two-fiuid model is given 

by (jAndersson. Glampedakis fc Samuelssonll2009t ) 



2 B 

va = 



47rpp 



(25) 



which depends on the proton density instead of the total 
mass density as in the single-fluid case. 

The perturbation of the gravitational potential is ob- 
tained by solving the linearised Poisson equation 



V5$ = A-kGSp . 



(26) 



Due to the symmetry of our background configurations, 
any perturbation variables may be Fourier expanded with 
respect to the coordinate (j). More precisely, in orthonormal 
spherical coordinates the mass density as well as any other 

J ierturbation func tion can be written in the following form 
Jones et al1l2002l ) 

5p— (t, r, 0) cos m<^ + (5pm (i, ''i ^) sin m<j!>] , (27) 

m=0 

where m is the azimuthal harmonic index. With this trans- 
formation the perturbation equations decouple with respect 
to m and the problem becomes two-dimensional in {r,6). 



2.3.1 Boundary conditions 

For non-rotating stars with purely poloidal or toroidal mag- 
netic fields, the two-dimensional numerical domain extends 
over the region Q < r < R and < < 7r/2, where R is the 
stellar radius. To study the evolution of non-axisymmetric 
oscillations (m 7^ 0) we must therefore specify the boundary 
conditions at the star's origin, rotation axis, equator and 
surface. In particular, for m > 2 the regularity of the lin- 
earised equations leads to a zero condition for all the scalar 
and vector perturbation fields at the origin (r = 0) and rota- 
tion axis (6 = Q), e.g. 5p = 0. At the equator (6 — 7r/2), the 
reflection symmetry spli ts the pertur bation variables in two 
sets of opposite parity (| Jones et al.l [20021. If we represent 
all the scalar and vector perturbation flelds by, respectively, 



the mass density and mass flux perturbations, the polar-led 
perturbations satisfy the following conditions at the equator: 

deSp = der = f = def^ = 0. 



The conditions for the axial-led class are 



sp = .r = dof = .r = Q 



(28) 



(29) 



The boundary conditions for the magnetic fleld pertur- 
bations depend o n the background fleld configuration 
iLander et al.|[20ldl ). For a poloidal field we have 



for polar-led perturbations and 

deb'' = fe" = deb'f' 







(30) 



(31) 



for the axial-led class. For a toroidal field, by contrast, the 
polar-led perturbations obey condition (|3ip and the axial-led 
perturbations obey pO[) . 

In our purely fluid models, i.e. with no crust, the two 
fluid components extend over the entire stellar volume. The 
crust is an important constituent of a neutron star which 
we intend to implement in a future work. In the current 
study the position of the star's surface (r = R) is given 
by the vanishing of the proton and neutron mass density, 
i.e. Pp — pn = 0. For the perturbation variables we require 
that the Lagrangian perturbation of the chemical potentials 
vanish at the surface, i.e. 



A^x = Sp^ + • V/ix = , 



(32) 



where ^x is the Lagr angian displacement a ssociated with 
each fluid constituent (|Andersson et al1l2004l '). These condi- 
tions may be expressed in terms of pressure and chemical 
potential perturbations as AP = A/3 = 0. For polytropic 
models, the surface condition m ay be further simplifled to 
SP = l|Passamonti et al.ll2009l '). From the deflnition of the 
stellar surface (px = 0) it follows that all the "flux" variables 
are zero at r = i?, i.e. f — D — b — 0. 



2.3.2 Numerical Code 

We study the time evolution of equations (|14|) - (|17| ) and (|23|) 
with a numerical code based on a MacCormack algorithm. 
The code is an ex tension of those a lr eady developed for a 
superfluid s tar by Passamonti et al] (|2009l ) and a magne- 
tised star bv lLander et al.h 2010l ). Iiithese two references the 
reader can flnd all the technical details. In order to speed up 
the numerical simulation we neglect the perturbation of the 
gravitational potential and adopt the so-called Cowling ap- 
proximation. In this way we do not have to solve the elliptic 
equation (|26p which is time consuming. This approximation 
has anyway a tiny effect on the Alfven mode frequencies and 
it is suitable for the aims of this work. 

We consider for the axial- and p olar-led pertur b ation s 
the same initial conditio ns used by iLander et al.l (|2010l ) 
an dLander fc Joned (|201ll '). 



3 RESULTS 

We are interested in studying the effects of two-fluid physics 
on magnetar dynamics. A realistic magnetar model has 
complex physics which includes superfluid/superconducting 
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Figure 1. Polar-led m = 2 Alfven modes for an unstratified 
star with a purely poloidal magnetic field with average magnitude 
B = 10^^ G. This figure shows the dependence of mode frequency 
on the proton fraction for a stellar model with zero entrainment. 
On the left vertical axis, the mode frequency u = a/{2iT) is shown 
in physical units for a star with M = IAMq and R = 10 km. 
The right vertical axis instead shows the dimensionless mode fre- 
quencies i//{Gpci)^^^ , where po is the central mass density of the 
background model. In the limit Xp I we recover the single-fluid 
results. The modes alread y determin e d (no t determined) in the 
single-fluid limit by iLander fc Jone i 1I2011I ') are denoted with a 
cross (flUed-circle). 



components, crust, realistic EoS, a strong magnetic field and 
magnetosphere. In this work we make a first step toward the 
introduction of superfiuid physics in magnetised stars, and 
study its impact on the Alfven mode spectrum. Here, we fo- 
cus on purely fluid models, while the presence of a crust and 
the study of magneto-elastic oscillations will be addressed in 
future work. 

We consider both stratified and unstratified stellar mod- 
els which can be obtained with an appropriate choice of the 
polytropic indices in the EoS (|10p . More precisely, models 
with constant (non-constant) proton fraction can be deter- 
mined by setting A'^p = A'n (A'p 7^ Nn)- For the unstratified 
case, we choose A'n ~ Np — 1 and study the effects of en- 
trainment and proton fraction on the oscillation mode fre- 
quencies. For the stratified case, we construct models with 
an increasing composition-gradient stratification and study 
its impact on the spectrum. 

This work studies non-axisymmetric oscillations of stars 
with purely poloidal and purely toroidal fields. We are able 
to investigate both axial- and polar-led modes of a toroidal- 
field star, but the axial-led perturbations of a poloi dal field 
are subject to the rapidly-growing Tayler instability l|Wrightl 
I1973I : iMarkev fc Tavlei|[l973 : iLander fc JonedlioTll ). which 
dominates the evolutions. In this case we consider only 
polar-led modes. We specialise to the case of oscillations 
of azimuthal index m = 2 for brevity. Whilst it is possible 
to study modes of higher m using our code, we anticipate 
that they will be similarly affected by two-fiuid physics; in 
addition, higher-m modes are probably more susceptible to 



Figure 2. Polar-led Alfven m = 2 modes for an unstratified 
star with a purely poloidal magnetic field with average magni- 
tude B = 10^" G. This figure displays the mode dependence on 
the entrainment e* for a star with a constant proton fraction 
= 0.1. A realistic parameter space for the core's entrainment 

1/2 

without strong pinning would be 1 < ej < 1-7 (see text). In the 
limit of large effective masses mj 2> nix, i-e. e* ~ Xp, the mode 
frequencie s tend as expect ed to the single-fluid results determined 
by Lander fc JonesI 1I2OIII) . The notation used in this figure is the 
same as in figure [T] 



damping in a real magnetar and hence less likely to survive 
to produce the observed long-lived QPOs. 



3.1 Oscillations of unstratified magnetars 



We begin by studying unstratified models, which are de- 
scribed by the EoS ^ with iVp = TVn = I. 

F rom the plane-wave analysis of lAndersson et al.] 
(j2009l '). we expect the Alfven mode frequencies to scale as 



(33) 



where ctq is the mode frequency of a single-fiuid star, and 
is given by 



1 - £n 



(34) 



1 - En 

The presence of the proton fraction Xp in equation (|33|) is 
due to the different definition of the Alfven velocity in a 
two-fiuid system, as shown in equation (|25|) . Equation (|33|) 
will guide the analysis of our numerical results, and we also 
expect that in the appropriate 'single-fiuid' limits (in terms 

in ILander et al, 



of Xp, £») we sh o uld re cover the mode frequencies reported 



iLander fc Jones 



2011 



( 201c ) for a background toroidal field and 



for a background poloidal field. 
For all the results reported here for unstratified stars, 
we take a fiducial average field strength of 10^^ G, which 
corresponds to a polar-cap value of 5 x 10^^ G for a poloidal 
field. This is a little larger than observed magnetar fields, 
but the correspondingly shorter Alfven timescale allows for 
faster numerical evolutions. We have, however, checked that 
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Figure 3. Axial- and polar-led m = 2 Alfven modes for 
an unstratified star with purely toroidal magnetic field with 
B = 10^^ G and zero entrainment. This figure shows the modes' 
sc aling with proton fr action. The mode frequencies determined 
bv lLander et alj l|201Cl ) are recovered in the single-fluid limit. See 
figure [T] for the notation used in this figure. 



for all our models the mode frequencies exhibit the expected 
linear scaling with magnetic field strength. 

In figure [T] we plot the dependence of polar-led Alfven 
mode frequencies on the proton fraction Sp, for a star with 
a purely poloidal field and zero entrainment. In this and fol- 
lowing figures, the mode frequencies v = o"/(27r) are shown 
both in dimensionless units, v / [GpoY^^ (where po is the 
central mass density), and in physical units, for which we 
have considered a star with mass M = IAMq and radius 
R = IQ km. It is evident from figure [1] that the expected 
relationship from the plane-wave analysis, i.e. a cc Xp , 
is borne out by our full numerical analysis. This is not sur- 
prising as in a non-stratified model the relevant parameters 
are constant. Most notably, we are able to determine the 
spectrum from a small proton fraction, Xp = 0.05, up to the 
single- fluid case Xp = 1, where we recover the results for a 
single-fluid magnetised star. 

Considering a model with Xp = 0.1 and a purely 
poloidal magnetic field, we show in figure [2] the scaling of 
the polar- led Alfven modes with the entrainment e*. We 
explore a wide parameter space: from the fully-entrained 
'single-fluid' limit, to entrainment values expected in the 
neutron star's core, passing by the unentrained case £* = 1. 

1/2 

The agreement with the expected scaling a oc ej is again 
very good. The single-fluid results can be obtained in the 
strong entrainment limit with large effective mass mj ^ m^, 
which leads to e* ~ Xp and thus to ct ~ (tq, see equation (|33p . 

The typical effective masses expected in a neutron star 
core without strong pinning range in the interval 0.93 < 
m-n/mn < 1 and 0.4 < mj/nip < 0.95 l|Chamell [20o3 v 

1/2 

These values, which lead to 1 ej 1.7, may increase 
the mode frequencies with respect to the single-fluid case 
by a factor of five in a star with Xp — 0.1. It is then 
natural to wonder whether this strong effect on the spec- 



trum may be influenced by the presence of an elastic crust. 
The conditions in the inner crust, where we expect a lat- 
tice of ions permeated by a gas of superfluid neutrons, 
can be in fact very different. Recent calculations show that 
superfluid neutrons may be efficiently entrained by nuclei 
due to Bragg scattering, and that the neutron effective 
mass at the bottom of the crust can be large as m* ~ 
14mn (|Chamel| |2012h. This strong entrainment can limit the 
relative two-fluid motion and produces a 10% correction of 
the shea r mode frequenc i es det er mined with a single-fluid 

; ISamuelsson & Andersson| 



mode l jA ndersso n et aU 200S 
iPass 



l2009l : iPass amonti fc AnderssonI [2012; Sotani ct al. 2012,). 



However, the crust's entrainment may be less relevant for the 
Alfven modes due to their global nature. A rough estimate 
can be determined by taki ng the effective masses of the crust 
and the core calculated bv lCarter et al.l (|2005l ') and lChamell 
(2005, 2006, 2008). If we insert th ese values in the unstrati- 
ficd model with Xp = 0.1 used by IPassamonti fc AnderssonI 
(2012) and take an average of the entrainment profile, we 
obtain ^/ (e*) — 1.2. In our purely two-fluid star this value 
would increase the Alfven mode frequencies by a factor of 
3.8 with respect to the single-fluid case. This estimate must 
be considered with caution, as the coupling of the magnetic 
field and the crust may lead to a different dynamical evo- 
lution and results. This issue will be addressed in a future 
paper. 

When the star has a purely toroidal magnetic field we 
are able to evolve stably both axial and polar initial data 
and extract mode frequencies. In figure[3]we show the effects 
of proton fraction on the axial- and polar-led Alfven mode 
frequencies for a star with a toroidal field and zero entrain- 
ment. In this case too the scaling of the mode frequencies 
is well approximated by equation H33p and we are able to 
recover the results of single-fluid magnetised models. 

From figures [T] and [3] we are able to compare the ef- 
fect of different magnetic field geometries on the oscilla- 
tion spectrum of a magnetar. Let us concentrate on the 



results for x. 



-1/2 



3.2, i.e. Xp 



0.1, and just the polar- 



led modes (ignoring the single axial-led mode in figure [3]). 
For a background poloidal field we find five widely-spaced 
Alfven modes in the broad range 50 — 500 Hz. By contrast, 
for a purely toroidal field there are only three modes, in 
the narrower interval 180 — 320 Hz. These differences sug- 
gest the possibility of constraining the field geometry with 
future QPO observations. 

Finally, in a star with a toroidal field we find the Alfven 
modes again scale in the expected manner with entrainment 
(cr oc ei ), as for the poloidal- field case. Since the plot of 
this contains no new information, however, we omit it for 
brevity. 

3.2 Oscillations of stratified magnetars 

We now turn our attention to the oscillation spectrum of 
stratified models and focus on purely toroidal magnetic 
fields. The poloidal configurations are currently less numer- 
ically stable, and we want to return to these at a later date. 
We use a weaker magnetic field than in the previous subsec- 
tion here, with an average strength of 5 x 10^^ G — again 
for reasons of numerical stability. 

The effects of stratification can be studied by using 
equation (|10|l with A'^n 7^ A'p. We consider a sequence of 
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stars with various composition gradients by setting Nn ~ 1 
and exploring a range of proton polytropic indices, 0.7 < 
A'p < 1.5. For this range we are able to run the code for a 
sufficient evolution time to extract mode frequencies reliably 
(15+ Alfven times). The central proton fraction is Xp — 0.15 
for all these models and tends to zero (unity) for A^p > 1 
{Np < 1). The frequencies of the axial- and polar-led m — 2 
Alfven modes are shown in figure |4] for stars with toroidal 
field B — 5 X 10^"' G and zero entrainment. The mode fre- 
quencies v — (t/(2-k) are given in physical units for a star 
with M = IAMq and R = 10 km. Note that in this strat- 
ified stellar sequence the central mass density po does not 
scale linearly with Np. Therefore, the dimensionless mode 
frequencies cannot be shown in figure |4] together with the 
physical values. 

Our results show that composition gradients may signif- 
icantly affect the Alfven modes with respect to unstratified 
models. The m = 2 modes at A^'p = 1.5 are a factor of about 
1.38 larg er than the unstratifi ed case A^p = 1. As shown re- 
cently bylLandereFaO (|2012l ). a " realistic" proton gradient 
can be reproduced by the EoS (|10|l with a proton-fluid index 
close to A^p = 2. Unfortunately, this value is slightly beyond 
our current numerical capacity. However, considering a lin- 
ear fit of the mode frequencies shown in figure |3] we find 
that a model with Ap = 2 should oscillate with frequencies 
roughly 1.67 larger than the Ap = 1 case. For instance, the 
linear fit for the polar-led Alfven mode which has the lowest 
frequencies in figure |4] is given by 



(20.0 + 56.9Ap) Hz. 



(35) 



From the unstra tified model we can identify the modes 
already calculated bv lLander et al.l (|2010l ) in the single-fluid 
limit. In figure |4] they are represented with a cross. 



4 DISCUSSION 

The main result of this work is that the multi-fluid 
physics of a magnetar is likely to result in considerably 
higher Alfven-mode frequencies than a standard single-fluid 
barotropic model would indicate. Although this was pre- 
dicte d at a qualitative level by an earlier plane-wave anal- 
ysis (jAndersson et all [20091) , this paper is the first quan- 
titative study of the oscillation spectrum of a multi-fluid 
magnetar. 

We have studied the time evolution of non- 
axisymmetric oscillations of purely two-fluid stars with su- 
perfluid neutrons and normal (not superconducting) pro- 
tons. The effects of an elastic crust will be included in fu- 
ture work. We have considered various magnetic field geome- 
tries, proton fractions, entrainment and composition strat- 
ifications and analysed their impact on the Alfven mode 
spectrum. 

Starting with unstratified models, we find that realis- 
tic values for the proton fraction and entrainment lead to 
Alfven mode frequencies three to four times larger than the 
single-fluid results. If we move from an unstratified star to 
one with a more realistic composition gradient the mode fre- 
quencies may be increased by another factor of ~ 1.67. From 
these results we can extrapolate that a 'typical' multi-fluid 
magnetar could have QPO frequencies which are roughly a 
factor of flve/six higher than expected from a single-fluid 
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Figure 4. Axial- and polar-led m = 2 Alfven modes for a se- 
ries of stratified stars with purely toroidal magnetic field with 
B = 5 X 10^^ G, zero entrainment and A'n = 1. The horizontal 
axis displays the proton polytropic index Np which is related to 
the degree of stratification. Stratification is zero for A^p = 1 and 
increases for values of A^p much different from TVn. The central 
proton fraction is Xp = 0.15 for all these stellar models. On the 
vertical axis we show the mode frequency in physical units for a 
star with M = 1.4Mq and R = 10 km. The curves are represented 
with the same notation used in figure [T] 



unstratified model. Based on our results and using equa- 
tions psp and p5[) . we may summarise the expected scaling 
of Alfven modes in the presence of superfluid physics as: 



6.3(70 



0.15 4-0.85 ( ^ 



Vl.3 



1/2 



fp(0) 
0.1 



-1/2 



(36) 

where a is a mode of a multi-fluid magnetar, ao is the corre- 
sponding single-fluid mode, and a;p(0) is the central proton 
fraction. 

Furthermore, our results suggest that there is some hope 
for constraining the magnetic fleld geometry of a magnetar 
from the QPO frequency distribution. For a typical unstrat- 
ified stellar model with a poloidal field we find five Alfven 
modes, in the wide range 50 — 500 Hz; instead, for a toroidal 
field there are three modes in a far narrower range, 180 — 320 
Hz. 

Let us consider a magnetar with average field strength 
B = lO^^G. Although this seems rather high, it corresponds 
to a more reasonable polar-cap value of 5 x lO'^^ G if we as- 
sume a poloidal-field geometry; if the field is more 'buried' 
(e.g. with a strong interior toroidal field) then it may have 
an average strength of 10^^ G but a polar-cap value of only 
10^^ G. For this field strength, typical single-fluid lowest- 
order Alfven mode frequencies are roug hl-fl 10 - 50 Hz for 
m = and 50 — 150 Hz for m = 2; hig her /-multipoles or 



2008 



values of m produce higher frequencies (iSotani et al.l 

IColaiuda et al] l2009l : iLander et al] l2010l : iLander fc Jones! 



This refers to m = axial modes; ISotani fc Kokkotaj ||2009| ') 
finds a fundamental polar m = mode at 300 Hz. 
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I2OIII : iGabler et al]|2012l ) . Temporarily ignoring the various 
simplifications in our model, let us look at the repercus- 
sions of formula p6p for identification of magnetar QPOs. 
The strong 150 Hz QPO of SGR 1806-20 or the 155 Hz 
QPO of SGR 1900-1-14 may be interpreted as an axisym- 
metric Alfven QPO of the magnetar's core — the multi- fluid 
equivalent of a 25-Hz peak predicted by single-fluid models. 
Similarly, the long-lived 625 Hz peak of SGR 1806-20 would 
then be the multi-fluid equivalent of a 100-Hz single-fluid 
mode, which is in the range expected for an m = 2 Alfven 
mode. 

The above discussion does not provide an interpretation 
for the lower-frequency magnetar QPOs and neglects the ef- 
fect of the crust. The effect of superfluidity on the crustal 
shear modes is similar to what is described in this paper for 
the A lfven modes. However, a realistic entrainment (jChamell 
|2012| ) may have a very different quantitative impact on these 
two classes of modes. In fact, the effective mass of neutrons 
may be quite large at the bottom of the crust leading to a 
total correction of the shear mode of about 10% with re- 
spect to the single-fluid models. The Alfven modes however 
may be less affected by this entrainment configuration which 
is confined only in a limited region of the star. In conclu- 
sion, multi-fiuid dynamics with a realistic entrainment may 
increase significantly the Alfven modes as described in the 
last paragraph, but only produce a 10 percent correction to 
the shear modes. Therefore, we suggest that the observed 
magnetar QPOs may fall into two classes: below roughly 
50 Hz they may be magneto-elastic crustal modes, whilst 
above this value they could represent Alfven oscillations of 
the multi-fluid core. 

The attractive feature of our interpretation is that it 
does not rely on the long-term excitation of high multipoles 
of the magnetar's oscillation modes, but instead suggests 
that observed QPOs are all low-order modes originating 
from the core or crust. Nonetheless, our description of a 
magnetar is still rather approximate. We must refine our 
superfiuid magnetised models and also include an elastic 
crust. One key issue is whether magnetically-modified elas- 
tic modes are really able to explain all the lower-frequency 
QPOs. In the axisymmetric case, two-fiuid physics may move 
the Alfven continuum to higher frequencies. In this way, 
more crustal modes could be outside the continuum and live 
for a longer time. Alternatively, some low-frequency QPOs 
could still be core modes if the magnetic field is weaker than 
expected, or if the two-fluid enhancement of axisymmet- 
ric oscillations is less signiflcant than predicted by equation 

dMl)- 

Our conclusions are based on a multi-fluid star where 
the dynamics of neutrons and protons is essentially decou- 
pled at linear perturbation order, interacting only through 
the entrainment. We assume that the magnetar's internal 
field strength is above the critical value at which supercon- 
ductivity is destroyed and hence that the protons form a nor- 
mal fluid, but this is far from certain. If instead the protons 
form a type-II superconductor, the oscillation spectrum will 
certainly be affected. The nature of the oscillations changes 
and their frequ ency is altered by a factor of ^ Hc\/ B, where 

^ 10^^ G (lMendelllll99i ). In addition, the presence of a 
magnetic force on the neutrons and coupling between proton 
fluxtubes and neutron vortices could lo wer the frequency of 
oscillations with respect to our results (|van Hoven fc LevinI 



I2OO8I : [ciampedakis et al.ll201ll ). These important issues de- 
serve more quantitative attention in future studies of mag- 
netar QPOs. 
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